Final-state effects in the radio frequency spectrum of strongly interacting fermions 
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We model the impact of final state interactions on the radio frequency spectrum of a strongly interacting 
two-component superfiuid Fermi gas. In addition to a broad asymmetric peak coming from the break-up of 
Cooper pairs we find that, for appropriate parameters, one can observe a sharp symmetric "bound-bound" 
spectral line coming from the conversion of Cooper pairs in one channels to pairs/molecules in another. 
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Radio frequency (RF) spectroscopy may become a pow- 
erful probe of the many-body state of a gas of cold atoms. 
One indicator of this potential is the sharpness of hyper- 
fine spectral lines (T): orders of magnitude smaller than the 
~100 kHz interaction strengths found in interacting clouds 
of lithium atoms. As cold atom experiments begin probing 
exotic states of matter, this separation of energy scales may 
allow one to detect subtle atomic correlations. 

Despite this optimistic viewpoint, there appears to be 
major theoretical holes in our understanding of RF spectra, 
even at a qualitative level. The primary difficulty is that 
when radio waves change the hyperfine spin of an atom, 
the entire many -body state needs to adjust. Including these 
final state interactions is a nontrivial many-body problem, 
akin to the one which must be solved to understand the X- 
ray spectra of metals [2|. These are hard problems - there 
is no generic prescription for including these final state in- 
teractions, rather the solution will depend on the system at 
hand. Here we present a variational calculation in which 
we calculate the radio frequency spectrum of a strongly in- 
teracting superfiuid two-component Fermi gas, including 
arbitrarily strong short range interactions in the final state. 
We find structure in the spectra which have not previously 
been theoretically described. 

We consider a model where a gas of neutral Fermions 
occupy two different hyperfine states |1) and |2). Radio 
waves drive a transition from |2) to a new hyperfine state 
1 3). Interactions between atoms in the three states are de- 
scribed by scattering lengths a i2 , ai 3 , and a 23 . 

Previous approaches to understanding the RF spectrum 
of a superfiuid Fermi gas have either neglected final state 
interactions |[3j|4j|5j|6ll21, or included them with sum rules 
HO |9j [TOll . diagramatics ATI , or energy arguments lfl2l . 
Many of these previous works focussed on trying to gain 
quantitative understanding of the spectra for very specific 
parameter values. Here we present a straightforward varia- 
tional approach which complements these other works: we 
focus on gaining a global picture of the qualitative struc- 
ture of the RF spectrum for all parameter values. While el- 
ementary, our formalism is quite powerful - as was kindly 
pointed out to us by Giancarlo Strinati, it is equivalent to 
the zero-temperature limit of the BCS-RPA theory used by 
Perali et al. ifTTl to explore the unitary limit. 



Final state interactions can qualitatively change the RF 
spectrum. For example, early experiments |[T3l with the 
three lowest hyperfine states of Lithium atoms [14] found 
that final state interactions nearly canceled out the contri- 
bution from interactions in the initial state. This is due to 
the proximity of three wide Feshbach resonances: B 12 = 
834G (AB 12 = 300G), B 13 = 690G (A5 13 = 122G) 
and B 23 = 811G (AB 23 = 2 2 2G). The scattering lengths, 
which diverge on resonance, all become large at the same 
time. When a\ 2 and di 3 are similar the interaction effects 
cancel [8|. Later experiments El QE [T71 Q3 [JJ) found 
nontrivial spectra - some of which were used as evidence 
of pairing in these gases. 

Here we use a BCS ansatz for the initial state, and a re- 
stricted set of final states which includes only "coherent" 
excitations: the quasihole created in the |2) state is re- 
stricted to have the same momentum as the excited atom. 
The neglected excitations will shift our spectral lines and 
introduce broadening, however our approach should cap- 
ture the qualitative features of the spectrum. Our approach 
is exact in the limiting case of a,\ 2 = a 13 , where the spec- 
trum is a single delta-function peak, unshifted from its free 
space value. Furthermore, when a 13 = 0, it reduces to 
the BCS result: a broad asymmetric peak corresponding to 
the breaking up of Cooper pairs. Generic ally we find a bi- 
modal spectrum containing both the delta-function and the 
broad continuum, though for some parameter ranges the 
delta-function lies in the continuum and becomes signifi- 
cantly broadened - disappearing as a 13 — > 0. These results 
mirror those found by Chin and Julienne in studying the RF 
spectrum of a single molecule ifTTll . Interpreting the delta 
function peak as a "bound-bound" transition, we encounter 
final state pairs in parameter ranges where no free- space 
molecules exist. 

We construct a "phase diagram" (figure [TJ in the final in- 
teraction (ai 3 fci?) versus initial interaction (o^^f) plane, 
delineating the regions where we expect to see a bound- 
bound transition and where we do not. We also plot the 
interaction values used in three experiments lfl6l [T9l l20l . 
The experimental observations are consistent with our pre- 
dictions; in lfl9ll no bound-bound transition is seen, while 
in ll20l one is always seen. 

We only consider the spectrum of the homogeneous gas 
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FIG. 1: Phase diagram showing for which values of initial (a\ 2 ) 
and final (013) scattering lengths the homogeneous RF spectrum 
contains a bound-bound peak. Scattering lengths are normalized 
via the fermi wavevector, kp — (3Tv 2 n) 1 ^ 3 , where n is the atomic 
density. Symbols: parameters from lowest temperature data in 
the experiments - references 1161 . II 1 91 and 1201 correspond to 
open diamonds, closed circles, and open circles respectively. 

(measurable tomographically |fT9ll ), and avoid discussing 
the trap averaged spectrum, where inhomogeneous broad- 
ening obscures the basic physics lfT2l . Furthermore, we 
restrict ourselves to the unpolarized case (ni = n 2 ) at 
T = 0. 

The role of the hyperfine energies ej (j=l,2,3), is most 
transparent if we make the Canonical transformation to the 
field operators ipj >k — > e~ iejt ipj^ k , to arrive at a Hamilto- 
nianft = H+H RF withi? = H 1 +H 2 +H 3 +H 12 +H 13 , 
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We neglect the H 23 term 112111 . The RF detuning is uj = 
Urf — (^3 — £2) an d the free particle dispersion is e k = 
k 2 /2m — fi, where m is the atomic mass, and we work 
in units where H = 1. The interaction parameters A^ are 
connected to the scattering lengths between states \i) 
and I j) by 

1 mV ^ m 

Ay At: a,, ^ k 2 

where V is the spatial volume. Of particular importance is 
that if A12 = Ai 3 , then [X, H] = 0, implying that the RF 
spectrum will consist of a single unshifted delta-function. 

We make the variational ansatz that the system is initially 
in the BCS state 

|GS) = J] (u k + v k 4^l_ k ) I vac) . (5) 

Minimizing (H — [iN), and neglecting the Hartree terms 
|p22l yields a ground state energy E GS = (H) = ^2 k {e k - 



E k ) + A 2 /A 12 , with v 2 



(E k - e k )/2E k , and u 



(E k + e k ) /2E k . The gap A obeys the equation 

1 



A12 



^ 2E k 

k k 



1 



(6) 



and the number of particles is N = 2^2 k \v k \ 2 . The 
excitations will be described by quasiparticle operators 
Ik = u k ip hk - v k ipl _ k and r/ fe = v k ip\ k + u k ip 2 -k with 
energies E k = y 7 \e k - fi) 2 + A 2 . 

Fermi's Golden rule states that given an RF detuning u, 
the rate of transitions / to a set of final states |/) is given 
by the imaginary part of a spin response function, 



I{u) oc ^2\(GS\H RF \f)\ 2 5(u; + E GS - E f ) 
f 

= Im (GS| X f =-X |GS) = hnR(u) 



(7) 
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where H = H — E GS . We approximate this response by 
restricting the final states to those containing a single quasi- 
particle. By definition this captures the "coherent" part of 
the response where most of the spectral weight lies 11231 , 
and aside from the BCS ansatz is exact when a 13 = a 12 or 
when a 13 = 0. Introducing the incomplete but orthornor- 
mal set of states 
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we replace the H in (|7]) by H = Yl k k 1 
where H kk , = (k\ H\k') = 5 kk , (E k + e k ) 
Matrix inversion gives 
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In terms of scaled variables u 
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E x (lo -E x -x 2 + l) 



As previously discussed, there are two contributions to 
the spectrum: a continuum from breaking up Cooper pairs, 
and a discrete peak from the conversion of a Cooper pair 
into a pair in the new channel. Mathematically, the con- 
tinuum comes from the poles of the integrands of the G's. 
Since we have severely restricted the available final states 
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FIG. 2: Theoretical radio frequency spectra at a\ 2 — > oo, where kp = {?>iT 2 n) 1 ^ and ep — k 2 /2m. The discrete peak on the right 
two figures correspond to a "bound-bound" transition between Cooper pairs in the various channels. Moving towards the left, that peak 
merges with the continuum. The height of the discrete peak has been scaled down for display; the actual spectral weight in the peak is 
shown in figure [4] 

and neglected Hartree terms, the location of this contin- 
uum is independent of the final-state interactions and corre- 
sponds to all ijj for which we can find a k with u = E k +€k- 
ie. u) > \/A 2 + fi 2 — /i. In a more sophisticated theory 
the threshold will depend on ai 3 . The discrete peak comes 
from the condition 
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This condition can only be satisfied for sufficiently strong 
interactions. In particular, if a 12 = oo, we find that this 
peak exists if ai 3 > or if ai 3 y / 2m/7 ;$ — 1 ■ In the former 
case the excited state can be thought of as a molecular pair, 
but in the latter case it is more analogous to a Cooper pair. 
When O12 = ai 3 , one can recognize that if one sets uj = 
in Eq. ([TTJ one recovers the gap equation. In this limit all 
spectral weight resides in this zero-detuning peak. 

In the BCS-BEC crossover, the chemical potential de- 
pends on interactions and density in a nontrivial way. Fol- 
lowing convention we will quote our results in terms of 
the Fermi momentum of an ideal Fermi gas with the same 
density as the one we are considering, kp = (S^n) 1 ! 3 . In 
figures[2]and[3](top) we illustrate the a i3 dependence of the 
spectrum by showing the spectral density in the case where 
ai2 = oo. As revealed by figures [T] [2] and [3] the pair- 
pair transition peak breaks off from the continuum around 
£Ji3&F = —0.7. As depicted by the solid line in figure [4| 
most of the spectral weight is in that peak close to unitar- 
ity. For < a^kp <C 1 (beyond the right hand edge of 
figures [3] (top) andQ the spectral weight shifts back to the 
continuum. In that regime, the wavefunction overlap be- 
tween the large |1)-|2) Cooper pairs and the small |1)-|3) 
molecules becomes negligible. Further, the spectral weight 
in the bound-bound peak goes to zero when it hits the con- 
tinuum (left edge of the curves in figure |4|, in agreement 
with IT7I . Figure [3] (bottom) is similar to figure [3] (top), 
except that we show the a V2 dependence of the spectrum 
keeping ai 3 = oo. 

We now discuss how the structure seen in figures [2] and 
[3] can be measured experimentally. First, since all physics 




FIG. 3: [Top] Spectral density as a function of final state scatter- 
ing length ai3 and RF frequency lu, when initial state interactions 
satisfy a\ 2 = oo. [Bottom] Same plot, except as a function of 
initial state scattering length a\ 2 and uj with ai 3 = oo. Spectral 
weight increases from white to black. 

only depends on the dimensionless parameters a^kp all 
interactions can be modified by changing the density. If 
a,\2 = oo then only the final state interaction is changed 
when the density is varied. This is a quite powerful knob: 
as seen in figure [T] two of the experiments |[T6l [T9ll have 
been performed very close to the point in Fig.[3](top) where 
a bound-bound transition emerges from the continuum. A 
modest increase in density should allow the observation 
of this feature. Alternatively, one can produce different 
final state interactions by changing which internal states 
one uses in the experiment ||20l . Restricting ourselves to 
Lithium |[T4l in states |1), |2) and |3), one has six com- 
binations: for example one could take a |1)|2) superfluid 
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Hazzard for various insights, and Kathy Levin for correct- 
ing our sum rule arguments. We found discussions with 
Henk Stoof valuable for thinking about diagrammatic ap- 
proaches to this problem. We are particularly indebted to 
Giancoli Strinati who sent us detailed notes showing how 
to connect this variational approach to the BCS-RPA. This 
material is based upon work supported by the National Sci- 
ence Foundation through grants PHY-0456261 and PHY- 
0758104. 



FIG. 4: Fraction of spectral weight in the bound-bound peak as a 
function of final state interaction strength 043, for different values 
of ai2- When the two interactions are equal, i.e., 043 = a\2, all 
the weight is in this peak. 

and excite |2) to |3), or one could take a |2)|3) superfluid 
and excite |2) to |1). Using the notation that ij — k repre- 
sents a \i)\j) superfluid where one excites \ j) to \k), tun- 
ing the ij interaction to resonance, and taking typical cen- 
tral densities |[T9l . typical final state interactions (1 / a^/ci?) 
are J24) 12-3:-0.94, 21-3:-0.19, 31-2:2.58, 13-2:2.14, 32- 
l:-0.86, 23-1:0.18. 

Given the focus on sum rules in recent theoretical studies 
ll8ll9l [T0l . we compute the zeroth and first moments of the 
resolvent R cxac t, defined in eq (|7]), and our approximation, 
-Rapprox, from eq Q . The \ k) basis is sufficiently large that 
the zeroth moment, or total spectral weight, is preserved by 
our approximation, 
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(GS| XlX |GS) = ^v 2 k 
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Y,\{k\X\G^)\ 2 = Y,vl 



(12) 



The higher moments of both -R CX act an d -Rapprox, on the 
other hand, are dominated by long tails which arise from 
the short distance structure of the BCS pairs. In the absence 
of final state interaction, for example, i? ap prox(w) oc lo~ 3 ^ 2 
for large u. With final state interaction, i^approxl^) oc 
u;~ 5 / 2 and our approach reproduces the first spectral mo- 
ment calculated in (8). 

Finally we comment on the quantitative accuracy of our 
approach. The largest approximation that we make is the 
neglect of Hartree-type interactions. These terms will shift 
the spectral lines. We also neglected incoherent processes 
which will broaden the spectra. Another important approx- 
imation was our neglect of finite temperature effects. These 
could be important in experiments where T/T F ~ 0.1. 
Despite the severity of these approximations, we believe 
that our calculation is valuable for its simplicity. It illus- 
trates the important physics in a transparent manner. 

We would like to thank Tin-Lun Ho, Mohit Randeria, 
and Wilhelm Zwerger for motivational comments, Wolf- 
gang Ketterle for discussions about the experiments, Kaden 
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